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We investigate the effect of a two dimensional electron gas on the linear optical properties of CdTe 
quantum wells. We evidence experimentally the high energy tail of the exciton and charged exciton 
resonances which depends on electron concentration. Based on that, we show that the scattering of 
electrons with excitons and charged excitons is needed to be included in the matrix transfer calcula- 
tions to describe the reflectivity spectra. We demonstrate by time-resolved reflectivity experiments 
the importance of electron distribution in the resonance lineshapes. 
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I. INTRODUCTION 

Past experimental studies have pictured the hnear 
optical properties of quantum wells in the presence of 
a moderate background electron density. They showed 
that the exciton resonance is strongly affected by the 
electron population. First of all electrons screen the ex- 
citons and reduce their oscillator strengtlP^. Then, they 
fill up the conduction band and contrib ute t o shift the ex- 
citon resonance towards higher energie^^Ell, Finally they 
collide with electrons and broaden their absorption lineal. 
Experimental comparison between exciton-electron and 
exciton-exciton scatterin g m echanism in bulk GaAd^, and 
in GaAs quantum wells^, showed that the exciton- 
electron scattering efficiency is one order of magnitude 
larger than exciton-exciton process. Both scattering pro- 
cesses are enhanced for the two-dimensional excitons as 
compared to bulk excitons. 

Linear optical spectrum of modulation-doped quan- 
tum wells also features a charged exciton (trion) reso- 
nance below the exciton resonanc^. Theoretical calcu- 
lations of the binding ene rgy of charged excitons were 
performed using variationaP^^HI] a nd ful l solution of the 
three-particle Schrodinger equatioiJi^^. They showed 
that only the singlet state is bound in the absence of 
magnetic field. Considering the lineshape of the trion 
resonance, it was argued that the momentum of the elec- 
tron initially present in semiconductor and used in the 
absorption process of a trion yields a low-energy tail on 
the trion resonanc d^^ l ^'^ l A simple derivation of the trion- 
electron scattering was also performed in GaAs quantum 
well^^, but the overlap of the exciton and trion reso- 
nances did not allow clearing identifying if such simple 
model leads to correct predictions. 

In this paper, we investigate the effect of a two- 
dimensional electron gas on the optical properties of 
CdTe quantum wells. Due to the large energy separa- 
tion between exciton and trion resonances, we are able 
to evidence experimentally the existence of a high-energy 
tail of both exciton and trion resonances. An electron- 
exciton and electron-trion scattering model, taking into 
account both elastic and inelastic scattering processes, 
show quantitative agreement with our experimental lin- 



ear spectra. These findings, even being yet far, can go 
in the direction of the modern theory of excitons within 
an electron gas, that argue that a trion is intrinsically 
a many-body object, made of a hole interacting with 
all electrons in the systenP^. We use also time- and 
spectrally- resolved pump-probe experiments to investi- 
gate the effect of scattering of electrons with trion and 
excitons on the nonlinear reflectivity spectra. 

The paper is structured in the following way: We 
give the sample characteristics in Sec. |IT] After pre- 
senting the measured reflectivity spectra of the sam- 
ple in Sec. |III| we show, in Sec. |IV[ that a calcula- 
tion without exciton- and trion-electron scatterings can- 
not possibly describe our system. Sec. |V] summarizes 
the exciton-electron and trion-electron calculation re- 
sults and we demonstrate that linear optical spectra of 
modulation-doped quantum wells are quantitatively de- 
scribed by a matrix transfer calculation including both 
exciton- and trion-electron scatterings. In Sec. VII using 



time-resolved pump-probe experiment we illustrate the 
effect of the electron population on the exciton and trion 
lineshapes. We end with the conclusion in Sec. 



VIII 



II. SAMPLE STRUCTURE 

We worked on a one-side modulation doped 
CdTe/Cdo.27Mgo.73Te heterostructure^^, containing a 
single CdTe quantum well of 8 nm between 50 nm thick- 
ness barriers. This structure was grown on top of a 4.5 
mm of CdTe which is deposited above a GaAs substrate. 
A remote donor layer of iodine with thikness of 49 Awas 
embedded in the cap layer 10 nm apart from the quantum 
well, yielding a population of about 4 x lO^'' cm~^ ex- 
cess electrons in the QW. The electrons from the donors 
can either fall in the quantum well or be trapped by sur- 
face states, both processes competing. In this regime, at 
temperature of 5 K, two strong resonances separated by 
3 meV arise below the band gap; they are attributed to 
heavy-hole excitons (1625.7 meV) and negatively charged 
excitons (trions, 1622.4 meV). A control of the electron 
density can be achieved by illuminating the sample with 
light more energetic than the energy gap of the barrier. 
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i. e. larger than 2.03 eV. In this case, we create electron- 
hole pairs in the barrier. Holes are attracted by elec- 
trons trapped in the surface states, while electrons are 
repelled towards the quantum well. By increasing the 
light intensity, we can increase the electron gas density. 
For instance, the practical electron concentration range 
from 4.3 x 10^° cm~^ to 1.2 x 10^^ cm~^ is attained. The 
thickness of the barriers was chosen in such a way that re- 
flectivity spectra are absorption-like and can be directly 
interpreted. 



III. 



CW REFLECTIVITY EXPERIMENTS 



We performed reflectivity measurements of the sam- 
ple at temperature of 5 K. A white light source was spec- 
trally filtered so that any high energy component (> 1.7 
eV) likely to modify the density of the electron gas in 
the quantum well was suppressed. The light was coUi- 
mated and focused on the sample. The reflection was 
then spectrally resolved by an imaging monochromator 
and recorded by a CCD camera. We kept the incident 
light intensity low and made sure we were working in 
linear regime. The light coming from a blue GaN light 
emitting diode (LED) was focused on the sample and 
completely covered the white light spot. By monitoring 
the intensity of the blue LED we could change the density 
of electrons in the well. The diffusion time of electrons 
from the barrier to the quantum well is macroscopic; it 
usually took 5 s to 15 s for the system to reach equilib- 
rium whenever the LED power was tuned. Reflectivity 
spectra were recorded for various electron densities. 
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FIG. 1: CW reflectivity spectra obtained at 5 K for dif- 
ferent electron densities ranging from 4.3 x 10^" cm~^ to 



different electron densities. For increasing electron den- 
sities, we clearly observe a blue sh ift of the exciton res- 
onance as reported in the literaturd^^^. It is due to the 
filling of the conduction band by electrons. Electron-hole 
pairs can only be photo-generated above the Fermi level. 
To a good approximation, the shift is linear in the Fermi 
energy Ep. Thus the exciton-trion energy difference is 
given hy E = Eq + aEp, where Eq is the trion binding 
energy and a is an empirical parameter. For CdTe quan- 
tum wells, a = l.O'P^. Assuming that Eq —1.7 meV, 
we estimate the electron concentrations in the quantum 
well for each spectrum in Fig. 1. The densities that we 
obtained range from 4.3 x 10^° cm~^ to 1.7 x 10^^ cm~^. 

We note that in the measured spectra, also a high 
energy tail of exciton and trion resonances develops as the 
electron density increases. We will show that this finding 
can only be explained by considering electron interactions 
with excitons and trions. 



IV. REFLECTIVITY CALCULATION OF THE 
SAMPLE 

As a starting point, we chose to use the linear sus- 
ceptibility proposed byl^liji the framework of density ma- 
trix theory. It includes both exciton and trion contribu- 
tions, has a rather simple form and takes into account the 
nonzero momentum of the excess electron in the quan- 
tum well used in the process of trion photogeneration. 
We emphasize that it neglects any scattering with elec- 
trons. It reads 



(1) 



where x'^ {^) ^^'^ X^i^) ^r^, respectively, the contribu- 
tions to the exciton and trion susceptibility. They are 
given by 



X- 



' uj - ujx - ilx' 



X^M =/t J dqrieiq) 



UJ — LOT + W^q — *7T 



(2) 
(3) 



where iox iy>x) is the exciton resonance energy, (/t) 
a contribution to the exciton (trion) oscillator strength, 
and 7x (7t) the exciton (trion) homogenous broadening 
mainly due to phonons; nj^q) is the Fermi-Dirac distri- 
bution of the electrons and Wq = h^q^Mx /'2,mf.MT is 
a correction to the trion resonance that comprises trion 
center of mass energy and electron initial momentum q, 
with TTie, Mx and the electron, exciton and trion 
effective mass, respectively. For a given initial electron 
momentum, the strength of the integrant is weighted by 
the optical matrix element 



Figure 1 shows the reflectivity spectra obtained for 



M'{q)^ / dp2V'u(0,P2)e' 



(4) 
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where rj = Mx/Mt and ipfsiPi^ P2) is the trion Is wave- 
function, with pi = Tie — Th the relative position of the 
first electron to the hole and p2 — T2e ^ ^'/i from the 
second electron to the hole. 

For the trion wave function, we relied on a simple 
variational function^^ with only two variational parame- 
ters A and A'. It is given by 



(5) 



Where A/r is the normalization factor. With this approx- 
imation, the light coupling element M{q) becomes 



Mriq) = 47Vt 



A7 



1 



A^(l + (ATfc)2)3/2 



\t ■*— > Ay 



(6) 



We evaluated, in Appendix C, the error induce by 
using the Chandrashekar's wavefunction. In Fig. [7] we 
compare the result obtained with variational functions 
|C1| and |C2[ They do not differ much and we conclude 
that Eg. |C3| is a reliable approximation. 

Exciton and trion resonance may also be inhomoge- 
nously broadened, due to quantum well interface rough- 
ness. Building on Andreani et al.^i^ work, we supposed 
that the broadening can be described by a Gaussian dis- 
tribution function. Thus, for excitons, we substituted the 
dielectric susceptibility Q by the convolution function 



1 



Av — exp 



inhom 



»7r/x|0x(O)P (u: 

inho™ ^ 



horn 



(7) 



is the complex error functior!^ and ^'^\^ara 

For trions 



the 



where w 

exciton inhomogenous broadening constant 
the same kind of equation holds, except the convolution 
is performed before integrating over q: 



inhom 



dqn,{q)\M^{q) 



LUt + Wq — ijT 
inhom 



(8) 
(9) 



Using the transfer matrix formalism described in Ap- 
pendix A and B, we calculated the mode of the field 
for our sample structure. The re fracti ve index and ab- 
sorption of the CdMgTe barrier^^IIIl CdTe quantum 
well, CdTe buffeil^SES GaAs substrate and iodine dopant 
were all implemented. For each resonance, we used the 
following fitting parameters: spectral position, oscilla- 
tor strength, homogenous linewidth and inhomogeneous 
linewidth. We fitted three densities: 4.9x 10^°, 6.7x 10^°, 
1.2x 10^^ cm~^. Results are shown in Fig.[2]and the value 
of the corresponding fitting parameters in Table |T] 

From the fits in Fig.|2] we see that it is impossible to 
reproduce both low energy and high energy tail of exciton 
and trion resonance because of their strong asymmetry. 




1620 1625 1630 1635 1640 
Energy [meV] 



FIG. 2: (red solid) cw reflectivity spectra measured at 5 K 
and for electron densities of 4.9 x 10^°, 6.7 x 10^° and 1.2 x 
10^^ cm~^. Best flt obtained with electron scattering (black 
solid) and without (blue dashed). 
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TABLE I: Fitting parameters of the reflectivity spectra for 
three electron densities tie. The homogenous broadening Fhom 
is essentially due to phonon interaction and corresponds to 
7x for excitons and 7t for trions (see Eqj2] an d [31 ). The 
inhomogenous broadering is defined in Eqs. [7[and [T] 



We note that in the measured spectra, the low-energy tail 
of the trion resonance predicted by'^ is marginal and can 
be fairly well reproduced by our fit. Conversely, a high 
energy tail of exciton and trion resonances develops as 
the electron density increases. Since we did not take into 
account the interaction between exciton and trion with 
electrons we cannot reproduce the correct exciton and 
trion lineshapes. In the next section, we will use a sim- 
ple model of exciton- and trion-electron interaction that 
will fit quite well the high-energy broadening of measured 
resonances. 
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FIG. 3: Top row: exciton-electron (a) elastic and (b) inelastic 
scattering. Bottom row: trion-electron (c) elastic and (d) 
inelastic scattering. 



V. NEUTRAL AND CHARGED EXCITON 
SCATTERING WITH ELECTRONS 

A theory of neutral and charged exciton electron 
scattering has been proposed by Ramon et al.'^. They 
investigated both elastic and inelastic scattering in GaAs 
quantum wells. They were able to fit their results on the 
high energy line shape of excitons, but failed to check 
their validity on the high energy tail of trions because in 
GaAs exciton and trion lines are not well separated. 

We applied their method to calculate the linewidth 
broadening due to the exciton-electron scattering in our 
CdTe quantum well. Here, we summarize the principles 
and show the results. 



the electron-scattering mechanism. This should be com- 
pared to an exciton linewidth of ~ 0.1 meV for acous- 
tic phonon scattering at T=5 K. It is explained by the 
exciton-electron interaction matrix elements that favor 
small energy-transfer transitions. At higher densities the 
effect of the phase-space filling becomes noticeable and 
effectively enlarges the exciton Bohr radius. Increasing 
the electron density further results in a shift of the max- 
imum linewidth from kx = to higher momenta. 

During the scattering process, the exciton can also 
ionized into a free electron-hole pair. The exciton 
linewidth rjj'^g*'*"^(fcj(:) due to this inelastic scattering 
is represented in Fig. |3jb). Although the magnitude of 
pineiastic^^^) is of the Same order as V^'^^^j^kx) , its func- 
tional dependence on the exciton in-plane momentum is 
very different. In particular, we note that the maximal 
linewidth is obtained at a very large momentum. This is 
due to the fact that in order for an exciton with initially 
small kx to be ionized, it must scatter on an electron 
with energy large enough to overcome its binding energy. 
At zero temperature this is only possible above the Fermi 
energy. 

The two scattering processes that we described 
above contribute to admix states with kx > to the 
kx ~ state. This admixture can be easily included in 
our previous calculation of the absorption by convoluting 
the imaginary part of the exciton susceptibility ([2| with 
a Lorcntzian function whose broadening is given by the 
sumrx-e(A:x) ^r^i-f (fcx) + ri^!!|''"=(A:x) (Fig.g- It 
becomes 

/ dLo'lm[xx{^-^')]^{uj\Tx^e{co')) (11) 

"'0 



A. Exciton-electron scattering 

An electron of momentum can interact elastically 
with an exciton kx- In that case, a momentum q is 
transferred from the electron to the exciton 



(fcjf , fee) — > {kx +q,ke- q) 



(10) 



In the process, the electron bound in the exciton can be 
exchanged with the free electron. The exchange scat- 
tering matrix element is very strong compared to di- 
rect Coulomb interaction. Using Fermi's golden rule, it 
is possible to calculate the scattering rate of the pro- 
cess 10 Summing over all final exciton states, results 
— in the first Born approximation — in the exciton 
linewidth T°^^^^^'^{kx) due to elastic scattering, as a func- 
tion of its initial momentum. 

Fig. ^a.) shows the computed values of T'^'^^^^^{kx) 
as a function of the exciton initial energy and the elec- 
tron density. The large linewidth obtained for rela- 
tively low electron density reflects the high efficiency of 



Use of the Kramers-Kronig relations yields the real part 
of the dielectric function. Since Tx-e{kx) is a decreasing 
function of kx, the Lorentzian peak is shifted to higher 
exciton energies. This is seen as the electron density 
increases. 



B. Trion-electron scattering 

For trions, the charge of the trion results in a diver- 
gence of its matrix elements in the limit of zero trans- 
ferred momentum. This divergence, originating from the 
infinite range of the Coulomb potential, is treated by ap- 
plying the Lindhard model for the potential screening. 
We use the Chandrashekar's variational function intro- 
duced above to perform the calculation. As for excitons, 
the broadening comes form direct and exchange scatter- 
ing. The broadening Tx--eikT) obtained is shown in 
Fig.^c) and|5];d) . 
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FIG. 4: A schematic picture of the convolution resuhing from 
electron scattering. The solid lines represent various exciton 
initial line shapes lm[xQW{iu — ^'))- The contributions of 
states with kx > to the kx ~ state are weighted by the 
value of the Lorentzian function ^x-el^^jT^x-eii-o)], given by 
the dashed line 
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TABLE II: Fitting parameters 



VI. ELECTRON SCATTERING CORRECTION 

Using the CdTe exciton-electron and trion-electron 
scattering calculation performed in the previous section, 
we were able to fit the three cw reflectivity spectra con- 
sidered in Fig.|2] For each resonance (exciton and trion), 
the free parameters were the spectral line position, the 
oscillator strength, the electron density, the homogenous 
and the inhomogeneous linewidth. The fits are shown in 
Fig. [2] and the parameters of the fit in Table [TTl 

The fits we obtained are exceptionally good. Only a 
very small inhomogenous broadening had to be included, 
which is consistent with the high quality of our quantum 
well. The homogenous broadening "fx and jt attributed 
to phonons remained quite small and more or less con- 
stant over the densities. This shows that we accounted 
for almost all electron induced broadening in our simple 
model. Had it been perfect, we would have expected the 
exciton and trion energy parameters to stay constant for 
all densities. Unfortunately, we could not get a perfect 
match of the spectrum lineshape by keeping the shift con- 
stant. Our model accounts for part of the shift but not 
all. We remind that it is a simple first order perturbation 
model and that we might expect higher order contribu- 
tions to play a role. The density obtained remained quite 
close to those calculated with a simpler model in Sec. |III| 

As for the oscillator strength, not much can be said, 
since we did not include any dependence on the electron 



density in our model. We did so intentionally because 
current theoretical models conclude that the trion oscil- 
lator strength depends on the volume of the sampld^. 
We think that it is a failure of the three-particle over- 
simplified trion model which is, strictly speaking, only 
correct in quantum dots. A correct description of the 
trion resonance should take into account the interaction 
of the excitons with all electron in the system. Such ap- 
proach seems to lead to trion oscillator strength that do 
not depend on the volum^^. 



VII. TIME-RESOLVED REFLECTIVITY 
EXPERIMENT 

Spectrally resolved pump and probe experiments 
were also performed in refiectivity geometry at 5 K in 
order to investigate the effect of the exciton-electron 
and trion-electron scattering on the nonlinear reflectivity 
spectra. We use the sample without illumination so with 
an excess electron population of about 4 x 10^° cm~^ 
in the quantum well. The 100 femtosecond output pulses 
from of a Ti:Saphire laser were split into two: a small por- 
tion circulated through a delay line and probed the exci- 
ton and trion resonances, while the major portion passed 
through a pulse shaper to generate 4 ps long and 0.5 meV 
wide pump pulses, with sufficient spectral resolution to 
selectively pump the trion resonance. The pump/probe 
intensity ratio was larger than 50 so that probe reflec- 
tivity spectra remained linear in the probe field. The 
polarization of both pump and probe pulses could be 
chosen independently to investigate the counter-circular 
CT+cr^ polarization configuration. We recorded the dif- 
ferential refiectivity spectrum AR — {R ~ Rq)/Rq for a 
given delay time between pump and probe pulses, i?o and 
R being the unexcited and excited refiectivity spectrum, 
respectively. 

We choose the excitation pulse to be at trion reso- 
nance, because by generating trions we induce inevitably 
a variation in the electron population. We decide to per- 
form the experiments with <T'^a~ circular polarization 
because with this configuration we avoid phase space fill- 
ing effects [34] enabling us to detect changing of the ex- 
citon as well the trion resonance lineshapes. In Figure 5, 
we show the differential reflectivity spectra obtained by 
pumping at trion resonance and probing with at delay 
time of 4 ps. At this short positive delay time, all co- 
herent processes have decayed awa}'^, being it then the 
best delay to observe nonlinearities on the exciton and 
trion resonances. Therefore, we can in this way check 
the changing in the exciton and trion resonances due to 
electron scattering. 

Indeed, the differential reflectivity spectrum in Fig. 
5 evidences variations in both trion and exciton reso- 
nances. We observe that the trion line is slightly red- 
shifted. We attribute this redshift to the reduced amount 
of free electrons in the quantum well by the trion gen- 
eration. This induces a decrease scattering contributions 
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from the electrons with the probed cr~ trions. Therefore, 
the trion high energy tail absorption will reduce and as 
a consequence, its lineshape will change shifting a little 
towards lower energies. The trion generation also brings 
changing in the exciton resonance. We observe the ex- 
pected bleaching of the high energy tail of the exciton 
absorption with an induced absorption of the exciton res- 
onance. With these results, by generating a trion popu- 
lation and therefore altering the electron distribution, we 
evidence the importance of electron scattering with exci- 
tons and trions in their lineshapes and in their nonlinear 
signal. 
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FIG. 5: Differential reflectivity spectra at delay time of 4 ps 
for a^a^ pump and probe configurations. The pump is 
at trion resonance. 



VIII. CONCLUSION 

In summary, we demonstrated that in modulation 
doped quantum wells, electrons play an important role 
on both exciton and trion resonance lineshapes. We ev- 
idenced experimentally that a high energy tail develops 
at both excitons and trions resonances as the excess elec- 
tron population is increased in the quantum well. We 
showed that with a model without considering exciton 
and trion-electron scatterings, we are not able to repro- 
duce the optical spectra. We used then a simple first 
order perturbation model taking into account elastic and 
inelastic scatterings of electrons with excitons and trions. 
The model is based on calculating the exciton - and trion- 
electron direct and exchange interaction matrix elements 
to derive the trion and the exciton lincwidth broadening. 
Finally, by including theses results in the matrix trans- 
fer calculations, we reproduced the highly asymmetrical 
resonance shapes of excitons and trions in the reflectiv- 
ity spectra. This is fundamental for the interpretation of 
the linear and non-linear properties of modulation-doped 
quantum wells. 
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APPENDIX A: MATRIX TRANSFER 
FORMALISM FOR A STRATIFIED MEDIUM 

The layered structure of our sample needs to be 
taken into account in order to consistently simulate 
the optical properties of the quantum well. We study 
the propagation of time-harmonic electromagnetic wave 
through a stratified medium comprising successive thin 
plane-parallel films. 

First we shall solve the wave equation 



V2£;(u;,r,,,z) 



^D{uj,r»,z) = 0. 



(Al) 



where the electric field E{u}, ry , z) and the electric dis- 
placement D{uj,r\\, z) are expressed in term of the fre- 
quency w, the in-plane position vector ry and the position 
along the growth axis z. 

For a homogenous dielectric medium, the elec- 
tric and displacement fields are linearly dependent 
D{uj,r\\,z) = e{ijj)E{Lo,r\\,z). Eq. (Al) takes the form 



of a Helmholtz wave equation for E{ijj, r\\ , z) 

2 , ^iy^y 



c2 eo 



E{uj,r\\,z) = 0. 



(A2) 



The dielectric dispersion may very well be complex. 
Because it is constant, the electromagnetic field is in- 
variant under in-plane translations (Bloch theorem) and 
solutions of ( A2 1 are plane- waves 



£;fcii(c^,r||,z) = efc„Ufc||(a;,z)e^'^i|-ni. (A3) 
Here fcii is the in-plane wave vector and the polar- 



ization vector. After substitution into (A2|, we are left 
with a one-dimensional problem for the mode function 



dz^ \ eo 



(A4) 



The solution represents two counter-propagating waves 
u(w,fc||,z) =^;^(fc||)e~*'=^"-H£;,.(fc||)e''=-^ (A5) 

|. (A6) 



/c^2 e(^) 



c2 eo 



The coefficient Ei and Ej. are two complex coefficient that 
have to be determined by imposing Maxwell boundary 
conditions at the interface between adjacent layers. Let 
us note that if e(w) is imaginary then kz is imaginary as 
well and the solutions are evanescent waves. 

We use transfer matrices to express the propagation 
of the electromagnetic wave at the interface between two 
consecutive films with refraction index ni and n2- Then, 
we combine the matrices to calculate analytically the 
mode of the whole stratified medium. 
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FIG. 6: Electric Fields at the boundary zi and z2 of a plane 
layer. 



In the framework of the one-dimensional prob- 
lem (A4) we define for each position z in space a two- 



dimensional vector containing the two coefficients in ( A5 ) 



Q 



El 



(A7) 



For a given structure, we would like to express the rela- 
tion between coefficient vectors at zi and Q*^^^ at 22- 
Because Maxwell equations are linear, Maxwell boundary 
conditions will result in a linear relation, which we write 
as 



Q(2) ^ MQ 



(1) 



(A8) 



The 2x2 complex matrix M thus defined is called the 
transfer or propagation matrix. If zi and Z2 were to 
be found in the same homogenous layer, from (A5) the 
transfer matrix would simply be 



horn 



gifej(2:2-2i) 

e- 







(A9) 



The Maxwell boundary conditions are invariant un- 
der time reversal. This means that the complex coeffi- 
cient of a transfer matrix M do not change if we reverse 
the time evolutior!^. Mathematically, we write this as 



M = TMT- 



(AlO) 



where T is the time reversal operator. We would like 
to calculate TQ. We remind that the amplitude of the 
electric field is given by the real part of its representation 
in terms of complex exponentials 



E{uj,k\\,r\\,z,t) ^-^[{Eie-'^'' 
The time reversal operator then acts as 

fE{uj, , r« ,z,t)^ E{u;, k« , r« , z, ~t) 



ifcii -rii 



3? [{E*i 



-ik^z _|_ ]7;*QikzZ-^ g-ife||-r||g-iwtj 



Thus, T reverse the sign of fcii —kn and acts on Q as 



T 







Ef 






e;_ 



(All) 



The time reversal invariance let us express the transfer 
matrix of a stratified medium in terms of its complex 
reflectivity and transfer coefficients. We consider the case 
of a unitary wave which is arriving on the left, reflected 
in the opposite direction with a wave amplitude r and 
transmitted to the right with an amplitude t. Then 



(A12) 



By applying the time reversal operator T on both side of 
the equality we get a second system of equations 



" t ' 




'Mil 


M12 




■ 1 " 







M21 


M22 




r 



' ■ 






M12 




r* 


t* 




M21 


M22 




1 



(A13) 



Noting that in terms of r and t, the reflectivity and 
transmissivity are R— \r\^ and T = a^2\^\^ with 





'3? 




ai2 = < 


3? 





(2) 



.(2) 



for TE polarization 
for TM polarization 



we solve the system of linear equations (A12l and (A13l 
and obtain 



(A14) 



M 



1 

t* 
r 



r 

t* 
1 

t 



(A15) 



Then, we use Eq. (Al5l here above to calculate the 



transfer matrix at an interface between two homogenous 
layer^^. We obtain the two transfer matrices for TE and 
TM polarization 



TE 
^12 



.(2) M 



k 



(2)- 



2k 



(1) 



2k 

.(2) 



(1) 



„(2) 



2k 



(1) 



2k 



(1) 



(A16) 



and 



M 



TM 



„TM 
"12 



,(2)fe(l)+e(l)fc(2) ,(2)^(1) _,(l)fc(2)- 



2VeWe(^ki^^ 2VeW^ki'^ 
e(2)fcW^eWfcf g(2)fcW+e(i)fci^^ 

(A17) 

These expression are quite general and are valid for com- 
plex dielectric constants as well. Combining those two 
matrices with the propagation matrix ( A9 ) we are able 



to calculate the electromagnetic modes in any multilay- 
ered structure. 



APPENDIX B: TRANSFER MATRIX OF A 
QUANTUM WELL 



We solve Eq. |Al| for a quantum well of with width L 
when TE polarized light propagates in the barrier. The 
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electric field is along the y-axis and ( Al I becomes 



rL/2 

{dl+kl)Ey+ex{^)p{z) / dz' p{z')Ey{z') - (Bl) 

J-L/2 



with k — \J e{uj) / tQU! / c. As p{z) is even for optical al- 
lowed transitions, it is convenient to treat separately solu- 
tions of the problem with definite parity under inversion 
of the z coordinate. 



For odd Ey symmetry, the integral in Bl is zero and 
there is no polarization contribution from the quantum 
well. The reflectivity takes the simple form 



(odd) _ _ik,L 



(B2) 



For even Ey symmetry, the solution of the second 
order differential equation (Bl) is 



Ey = A{uj,k;)cos{k^z)+ dz' G{z,z')p{z'). (B3) 



In this expression, the first term is the general even solu- 
tion of the problem {d1 + k1)Ey — 0, with A{uj, kz) a con- 
stant to be determined. The second term is a particular 
solution of ( B 1 1 and was simply built from the definition 
of the Green function {d^ + kl)G{z, z') = S{z — z'). We 
choose 

G{z,z') = -^sm{kz\z-z'\) (B4) 

which has the advantage to yield an even solution of Ey. 
Replacing ( B3 1 into (|B1[) gives 



(B5) 
(B6) 



rL/2 

Qi^z) = / dz p{z) cos{kzz) 

J-L/2 

.L/2 .L/2 ^ 

/ dzdz' —— sm{kz\z — z'\)p{z)p{z') 
'L/2J-L/2 2^2, 

(B7) 



Pikz) 



We consider Maxwell boundary condition at the interface 
z — —L/2. On the left side of the quantum well, we 
know that there are two counter-propagating plane waves 
propagating in the barrier given by Eqs. [X|and [X5l The 
Maxwell boundary conditions give 



V ' y 



L/2- 



UE'^'+4'')\z = -L/2- 

^ikz{E^^ -E^;^) 



E,, 



y\z = -L/2+ 



V '\z = -L/2- ~ '^''■^v\z=-L/2+ 



(B8) 



Substituing Eq. B3 in those two equation, we obtain the 



reflectivity for the even solution: 



(even) -^y 



E„ 



^dEy 



' TE 



ikzL 



E,, 



y ' 
A{(jj,kz 



ih^^Ey z = -L/2 + 

- iQ{kz)/2kz 



A{uj,kz)+tQ{kz)/2kz 
„tk,L 1 - X{^, kz){ia{kz) + Pju;, kz)) 
1 + X{'^,kz){ta{kz)- P{uj,kz)) 



(B9) 



where a{kz) — Q{kz)'^ /2kz. The total reflectivity is given 
by the mean of the odd |B2| and even reflectivity |B9| 



' TE 



1 /„(cvcn) (odd)\ 

2 \^'TE "1" 'TE J 

-ia{kz)x{^-,kz)e''''^ 
I + x{oJ.kz){ia{kz) - P{uJ,kz)) 



(BIO) 



The same kind of calculation apply for the transmissivity, 
which reads 



-i _ i / .(even) 
''TE ~ 2 I ''TE 



(odd) 
TE 



{\-P{u:,kz)x{^,kz))e 
1+ [ia{kz) ~ P{uj,kz)]x{uJ,kz) 



(Bll) 



If we now suppose that the quantum well has an infinite 
potential, the calculation of a and P functions becomes 
trivial. We obtain a ^ 1 and P ^ 0. We are left with 
the simple reflectivity and transitivity 



' TE 



ik^L 



-ix{uj,kz)/2kz 
l + ix{uj,kz)/2kz 



™ l + tx{io,kz)/2kz 



= 1 



' TE 



(B12) 



(B13) 



When used in Eq. |A15| these two expressions give the 
transfer matrix of a quantum well. The only information 
that is required to simulate the linear optical spectrum 
of our quantum well is the susceptibility. 



APPENDIX C: VARIATIONAL 
WAVEFUNCTIONS 

For the trion wave function, we relied on a simpler 
variational function^^, given by 



V't(Pi,P2) = {e 



X (1 + cR) 



-sR 



l + d{R- Ro) 



(CI) 



where R — — r2e- The variational parameters a, b, c 
d, s, Rq were obtained by fitting the numerical solutio n of 
the Schrodinger equatioiJ^. They are given in Table III 
for both positive and negative trion in units of the bulk 
CdTe bulk Bohr radius (7.75 nm). The overlap with the 
numerical wave function is excellent (0.9981 for X~). 
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a h 


c 


d 


s 


Ro 


X' 


1.1723 0.5471 


0.4809 


0.0682 


0.0013 


1.1614 


x+ 


1.1774 0.4294 


2.3943 


0.08431 


0.26173 


0.7393 



suit obtained with variational functions CI an d|C2[ They 
do not differ much and we conclude that Eq. |C3| is a re- 
liable approximation. 



TABLE III: Variational parameters for Sergeev and Suris' 
trion variational function 

Even if this variational function is analytical, its 
Fourier transform is not. To speed up calculations, we 
resorted to a simpler Chandrashekar's wave functioiJ^ 
with only two variational parameters A = 6.8 nm and 
A' = 15 nm. It is given by 

V't(pi, P2) = ^^T (eP^/^'P-/^' + eP^'^'-P^/^) , (C2) 

Where A/r is the normalization factor. With this approx- 
imation, the light coupling element M{q) becomes 

(C3) 

We evaluated the error induce by using the Chan- 
drashekar's wavefunction. In Fig. [7] we compare the re- 





Chandrashekar's 




Stebe's 











0.5 1 1.5 2 2.5 3 



q/Xx 



FIG. 7: Trion optical coupling as a function of the momentum 
q of the initial electron. 
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